Stochastic contribution to the growth factor in the ACDM model 
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We study the effect of noise on the evolution of the growth factor of density perturbations in 
| the context of the ACDM model. Stochasticity is introduced as a Wiener process amplified by 

an intensity parameter a. By comparing the evolution of deterministic and stochastic cases for 
different values of a we estimate the intensity level necessary to make noise relevant for cosmological 
tests based on large-scale structure data. Our results indicate that the presence of random forces 
underlying the fluid description can lead to significant deviations from the non-stochastic solution 
at late times for a > 10~ 3 . 
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O ■ I. INTRODUCTION 

Oh' The problem of describing the growth of small density perturbations in the universe consists of taking the differential 
equations of the fluid dynamics, rewriting them as a unique wave equation for the density contrast, and finding the 
solutions of this equation for different cosmologies. This is a classical problem of theoretical astrophysics, which 
began with the work of Jeans [1] , its general realtivistic generalization was performed by Lifshitz [2] . For a modern 
presentation of this problem see for instance [3]. According to the traditional approach, after the generation of the 
initial spectrum of density perturbations in the early universe, the subsequent evolution of clustering is deterministic, 
and does not admit a noise term in the dynamical equations. However, stochasticity can be present due to processes 
hidden by the coarse-grained description of the fluid and whose typical time and length scales are much shorter 
than those considered for large scale structure formation (see e.g. [4]). Actually, noise can be relevant in non- 
exceptional conditions, implying that even if it is very weak in the beginning, its effects can be amplified and have 
a non-negligible impact on the later dynamical evolution of density perturbations [5,6]. At the same time, a very 
basic aspect of stochastic phenomena is related to the differential law itself. It is well known that in a wide class 
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of real dynamical situations, the relationship between random excitations and the response of a physical system is 
governed by differential equations which cannot be strictly deterministic, but are stochastic (at a certain level) in 
nature (see e.g. [7]). This suggests that exploring stochastic effects on the evolution of density perturbations can be 
very important to modern cosmology, since the growth of large-scale structures is supposed to be used to constrain 
the dark energy equation of state, see e.g. [8,9]. Therefore it seems important to understand if the presence of a 
stochastic term in the differential equation still allows us to contrain cosmological models via large-scale structure 
data. In this work, we study the effect of noise on the evolution of the growth factor in the context of the ACDM 
model. Our aim is just to compare the evolution of deterministic and stochastic cases and estimate if the typical 
differences are similar or not to differences expected from assuming different cosmological models. 

II. GRAVITATIONAL EVOLUTION OF STOCHASTIC PERTURBATIONS 

We consider a generic instability scenario in which perturbations are generated by some mechanism in the early 
stages of the universe evolution, and start to grow under gravity when non-relativistic matter begins to dominate the 
density of the universe. The fluid is supposed to be pressureless and ideal, where particles do not cross. The equations 
describing the motion of the fluid are the continuity, Euler and Poisson equations. In Eulerian formulation they are 
respectively 

dp a 1 . 

Ji + 3 p+ v H = o, i 

at a a 
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dv 1 , a 1 

— + - v-V v + -v = V0, 2) 

at a a a 

V 2 (f> = 4nGa 2 (p- p b ), (3) 

where p?, is the homogeneous (background) density of the fluid and all the quantities are functions of the comoving 
coordinates x = r/aif), with a(t) given by the solution of the Friedmann equation, see for instance [3]. If we define 
the density field as 

p(x,t)=pi(t)[l + <f(x,t)], (4) 
where 5(pc,t) is the density constrast, the fluid equations become 

§ + V(l + *)v = 0, (5) 
<9v 1 , _. d 1 _ , 

— + - v-V v + -v = — V0, 6 

ot a a a 

V 2 = 4nGa 2 S. (7) 
Combining these equations, we find the second order differential equation for S 

^ + 2~-4nGp b 5 = A, (8) 
at 1 a ot 

where 

In the linear approximation, the term A is not important and we simply have 

2*^-4^ = 0, (10) 
ot z a Ot 

(see for instance [10]). In order to solve this equation, the following decomposition is useful: 

<y(x,i) = <y,(x)D(t). (ii) 

Since (|10p does not depend on spatial derivatives, S evolves only in amplitude, preserving its original shape in the 
linear regime. At the same time, the function D(t) satisfies the following equation 

D + 2-D - 4:nGp b D = 0, (12) 
a 

whose solution depends on the specific cosmological model adopted. This is the standard way to find how small 
perturbations grow in the post-recombination expanding universe. The fluctuations <5(x, t) clearly build a random 
field in space, since at a specific position x$ we cannot know the exact value of S. In most models <5 s (x) is supposed 
to have a Gaussian distribution, while D(t) is always a deterministic function evolving with time according to (|12| . 
Actually, nothing is asserted about D(t), which means that there is an implicit assumption in (TTTJ) saying that D(t) 
does not have a stochastic nature. In this work, we assume instead that noise can be present in the fluid due to the 
graininess of the underlying physical system of particles [11]. For a wider discussion on the possible origins of random 
forces in the large scale strucuture formation see e.g. [4,5,6]. A simple way to introduce a stochastic term in (TT2| is 
to suppose that the field is submitted to a zero-mean randomly fluctuating frequency C(*)j such that 

D + 2-D - [4irGp b + C(t)]D = 0, (13) 
a 

where £(£) = aW is the magnification (by a factor a) of a Wiener process (W), a quantity which is supposed here to 
embrace all the possible sources of random effects acting on the dynamical equation (TT2")) . Because of the stochastic 
nature of (p~3|) . we are interested in the statistical distribution of its solutions. This can be done directly by solving 
the equation numerically many times with the same initial condition and averaging over the results. The formal 
development of the theory of stochastic differential equations can be found in [12], while applications to stochastic 
phenomena in astrophysics are discussed in [13]. 
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III. STOCHASTIC CALCULUS 

In this section, we describe the method used to solve numerically the stochastic differential equation (SDE) (| 13|) . 
In short, we transform this equation into a set of first-order equations, take their integrals and iterate them from ti 
to ti + h, then from there to + 2h, etc. For sufficiently small values of h, this approaches the atual integral. The 
difference here is that instead of having usual integrals, we must deal with Ito (or Stratonovich) integrals, as described 
in [14,15]. Suppose that we have to solve the following SDE: 

d5(t) = m[t, 5(t)} dt + a[t, 5(t)]£(t) dt (14) 

where S(t = U) = Si and £(t) is the stochastic term. Here we assume that £ is the derivative of a Wiener process 
W(t). This is a continuous-time random walk with random jumps at every point in time, i.e., a step by step process 
in which the succession of steps is random and that in every step (jump) the variable changes values in a stochastic 
(random) way. The Wiener process W(t) is characterized by three facts: (i) W(0) = 0; (ii) W(t) is almost surely 
continuous; (iii) W(t) has independent increments for < s < t, e.g. [14]. Now, Equation (14) can then be re-written 
as 

dS{t) = m[t, 5{t)} dt + a[t, 5{t)} d[W(t)]. (15) 

(This is far from being a trivial change, as the Wiener process may not admit a time-derivative, but we ignore it here 
and follow the very basic approach in the numerical procedure) . Integrating both sides of the equation, we obtain an 
ltd integral, 

S(t) = Si + jf m(s, S(s)) ds + ^0, Hs)) d[W{s)} (16) 
whose iterative counterpart in the interval [tj,tj*_|_i] is 

/■tj+i 

S(tj + i) = Sj + / m(s,S(s))ds 
Jtj 

+ a(s,S(s))d[W(s)}. (17) 

Once in the form of an Ito (or rather, for technical reasons, a Stratonovich) integral, the equation is expanded about 
ti and solved with a fixed number of terms in the expansion (here, we included terms of up to three nested integrals), 
see chapter 5 of [14] for details. 

IV. GROWTH RATES AND COSMOLOGY 

Observational evidence suggests that the universe is currently accelerating, which would imply the existence of a 
significant unknown (dark) energy component pervading the whole universe as a homogeneous fluid [16,17,18,19]. To 
the moment, the existence of a cosmological constant seems to be the most economical explanation for the present 
data, although we have no fundamental physics solution to the coincidence and fine tuning problems, e.g. [20]. 
Trying to understand what is really behind the universe acceleration, a diverse set of projects designed to probe 
the nature of dark energy are in progress now. In fact, the increasing quantity and quality of large-scale structure 
data will soon allow us to use structure formation to constrain the dark energy equation of state, e.g. [21]. Dark 
energy affects the growing mode of density perturbations through the damping term H = a/ a in Eq. (13), making 
structure formation an interesting cosmological discriminator. However, differences in linear growth rates (in dark 
energy models) compared to a constant A for z < 3 are so small (of order 10 -2 ) that it may be extremely difficulty to 
rule them out observationally [22] . In this work we study the effect of a hypothetical stochastic contribution to the 
linear growth factor. Assuming it as a Wiener process, we would like to know the intensity level of this noise term 
to produce significant differences between stochastic and non-stochastic ACDM model. Accordingly, we take just the 
case of a flat universe with cosmological constant, so that 



(18) 
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a 




a f,ns 


1 A 

|A a , ns | 


0.0001 
0.001 
0.005 
0.01 


1.070 x 10 _b 
3.230 x 10~ 3 
6.040 x 10 -2 
3.200 x 10+ 2 


1.777 x 10" b 
3.287 x 10~ 3 
6.611 x 10~ 2 
3.233 x 10+ 2 


2.347 x 10" b 
5.830 x 10" 2 
3.584 x 10" 1 
0.668 x 10 +1 



TABLE I: Statistics for 200 functions in the entire time range. 



Using our simulation code, we performed a number of simulations for the case VL m = 0.3, Q\ = 0.7 and h = 0.7. 
We carry out 200 simulations for each of four values of a: 0.0001, 0.001, 0.005 and 0.01, within the redshift range 
z = 1000 — ► or t = — > 860 in our simulations units. The range of the magnification factor a is arbitrary, but 
chosen to be large enough to probe the levels at which the noise term turns important to the growth rates. Also, our 
choices of a respect the fact that we do not expect large signatures of the stochastic term near the last scattering 
surface due to the high entropy per baryon by this epoch [23]. The values of a studied here just produce late time 
(significant) effects. The quantity that we are interested in this analysis is the fractional error with respect to the 
growth rate of the ACDM model 



t _ A-D _ D stoc h — Dacdm , 1( ,n 

/ = -77- p; y w ) 

L> UaCDM 

Table U shows the effect of increasing a. It presents three basic quantities: cr/ ja is the largest standard deviation 
of all functions when compared with the average function; it indicates how much the functions spread around the 
average. <xy jns is the same, but compared to the non-stochastic solution rather than the average. Finally, |A ajns | is 
the absolute value of the largest difference (for all t) between the average of 200 simulations and the non-stochastic 
solution. As expected, the maxima happen at z = for all cases. In Figure 1, we present the evolution of AD/D 
for each value of a. Note that in all cases the presence of the stochastic term is not important at early times, but it 
becomes increasingly relevant at late times. We also see that the intensity level of the Wiener process can drive the 
growth rate contrast into a very noisy regime, with amplitudes higher than 10 -2 for a > 10 -3 . 



V. DISCUSSION 



Our stochastic simulations suggest that the presence of random forces under the cosmic fluid description can lead to 
significant deviations from the non-stochastic solution for the growing mode at late times. In fact, for intensity levels of 
a Wiener process higher than 10 -3 , the stochastic contribution to the growth factor would be of the same order as the 
predictions of typical differences between ACDM model and dark energy models [22]. Hence, this noise contribution 
would imply an additional difficulty to rule out models using large-scale structure data. In the present work, the 
introduction of stochasticity is completely ad hoc and our toy model based on the amplification of a Wiener process 
is quite arbitrary. However, our purpose is just to figure out if such contribution could be relevant to cosmology even 
for a simple implementation of the noise term. Actually, if noise should enter or not in the cosmic density evolution 
description is still a matter of debate, see e.g. [4]. Our results just raise the problem of distinguishing cosmological 
models using structure formation in the context of a cosmic fluid with underlying random forces. Of course, to pursue 
this effect more properly, theoretical efforts are necessary in order to understand the expected noise levels ab initio 
and the correct way to model and introduce them in the fluid description. 
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FIG. 1: AD/D against t (simulation units) for a = 0.0001, 0.001, 0.005, 0.01 (from top to bottom, following pannels (a) to (d)). 
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